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ABSTRACT 

We present predictions for the clustering of Extremely Red Objects (EROs) in a A cold 
dark matter Universe, using a semi-analytical galaxy formation model in combination 
with a cosmological N-body simulation. EROs are red, massive galaxies observed at 
0-''' ^ ^5 3, and their numbers and properties have posed a challenge to hierarchical 
galaxy formation models. We analyse the halo occupation distribution and two-point 
correlation function of EROs, exploring how these quantities change with apparent 
magnitude, colour cut and redshift. Our model predicts a halo occupation distribution 
that is significantly different from that typically assumed. This is due to the inclusion of 
AGN feedback, which changes the slope and scatter of the luminosity-host halo mass 
relation above the mass where AGN feedback first becomes important. We predict 
that, on average, dark matter haloes with masses above lQ^^h~'^ Mq host at least one 
ERO at 1.5 < z < 2.5. Taking into account sample variance in observational estimates, 
the predicted angular clustering for EROs with either (R — K) > 5 or (z — K) > 4.5 
is in reasonable agreement with current observations. 



1 INTRODUCTION 



Hierarchical models of galaxy formation have enjoyed many 
successes at both low and high redshifts. However, such 
models are still challenged by the large number of galax- 
ies observed with old stellar populations, many of which 
were already in place at redshifts higher than 2 = 1 
when the Universe was less than half its present age (e.g. 
ICimatti et"al] l2002al : IVaisanen fc JohanssonI 120041 ') . Obser- 
vationally it is possible to select such galaxies by requir- 
ing them to have very red optical and near-infrared colours. 
Extremely Red Objects (EROs) are commonly classified as 
galaxies with co lours redder than (R — /(')vcga ~ 5 (see e.g. 
iMcCarthvl 120041 ') . This classification was origi nally thought 
to se lect galaxies at very high redshift (z > 5, lElston et al.l 
Il988l ). However, later spectroscopic confirmation showed 
that E R Os were actually at 0.8 S, z < 2 (Cimatti et al.l 
l2002bl l3: IConseUce et 111120071 : IWilson et al1l2007l ). 

Currently there are no direct measurements of the spa- 
tial two-point correlation function of EROs, due to the 
small size of available surveys. The constraints on the clus- 
tering strength of EROs come from deprojecting angular 
clustering. Assumptions are required to extract the spatial 
correlation length from such measurements. Nevertheless, 
a consensus has been reached in that the observed clus- 
tering of EROs appears to be strong, with inferred cor- 
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these objects are more strongly clustered than the under- 



lying dark matter at the typical median redshift of these 
samples, which is close to z = 1. Improved measurements 
of ERO clustering require larger areas of the sky to be 
mapped. This is being achieved by various near-infrared sur- 
veys that are pushing ERO clustering measurements ou t to- 
wards larger angular pair separations (|Kim et al.ll201lh . 

Measurements of the clustering strength of EROs and 
its dependence on properties such as luminosity and colour 
will allow us to pin down the masses of the dark matter 
haloes which host these galaxies and hence to constrain 
the physics behind their formation. The high clustering 
amplitude for EROs implied by current observations sug- 
gests that these galaxies reside in massive haloes, in agree- 
ment with the predictions in the first paper in this series 
(I Gonzalez- Perez et al. I l2009l . hereafter Paper I). The cluster- 
ing signal also provides new constraints on the form of the 
feedback from both supernovae and active galactic nuclei 
(AGN) beyond those provided by the abundance of EROs. 

The clustering of massive red galaxies has pre- 
viously been interpr et ed u sing empirical models. 
iMoustakas fc Somerviuj (|2002l ) used the halo mass 
function and a power law form for the halo occupation 
distribution (HOD) to generate a simple model tuned 
to reproduce o bservations of EROs with H < 20.5 and 
{I~H)>3 bv lFirth et al] lj2002h . [Tinker et al.l (|2010l l also 
used a halo occupation distribution model to interpret the 
clustering of distant red galaxies (DRG), which typically 
have a higher redshift {z ~ 2) than EROs. Such models 
allow one to estimate the halo mass in which red galaxies 
are found, without addressing the physics behind how 
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the galaxies formed or why they have red colours. These 
calculations depend on assumptions about the form of 
the halo occupation distribution. As we shall show, this 
is a difficult proposition as the HOD of EROs is different 
from that of less ex treme, better studied galaxies (see also 
lAlmeida et"al]|2010l ). 

The abundance of EROs has also been tack- 
led in more ambitious calculations which do attempt 
to address the origin of EROs a nd thei r pho t omet- 



ric properties (iNagamine et aL 12005 : Kan g et al, l2006l : 
iFontanot fc Monaco! I2OI0I ). iNagainine et all (|2005l ) used 
mesh and particle based gas dynamic simulations to study 
massive galaxies at the redshifts of EROs. The simulation 
boxes used were relatively small, as demanded by the need to 
attain reasonable spatial and mass resolution in the calcula- 
tions. This makes it difficult to study rare populations with 
strong clustering using this ap proach. In order t o match the 
observed abundance of EROs, INagamine et al.l were forced 
to apply, by hand, quite a high level of extinction to all of 
their galaxies. Thus, all of their EROs are dusty star form- 
ing galaxies by construction, with no information about the 
number of galaxies with red colours due to their having old 
stellar populations. Semi-analytical galaxy formation mod- 
els have also been applied to study EROs. A key consider- 
ation here is that the model should reproduce the proper- 
ties of galaxies in the local universe, which are us ed to set 
most or all of the model parameters (|Baughll200^ ). Of par- 
ticular relevance for ERO predictions is the present day K- 
band luminosity function. The bright end of the luminosity 
function is dominated by early-type galaxies (jNorberg et al.l 
|2002| ). The challenge is for a model to match the abun- 
dance of EROs at intermediate redshift whilst still produc- 
ing the observed n umbe r of bright galajcies locally. Both the 
iKang et al] (|2006l ') and iFontanot fc Monaco! l|2010l ) models 
overpredict the number of bright galaxies in the K-band 
today. In particular, the IFontanot fc Monaco! l|2010l ) model 
overpredicts the number of galaxies by a factor of 10 at 
Mk — 51og h — —25.8, which corresponds to the knee of the 
luminosity function. 

In this paper we make predictions for the clustering 
of EROs ba s ed on the published semi-analytical model of 
iBower et al. I (I2OO6') in its original form. This model repro- 
duces many observations, including the present-day K-band 
luminosity function and the number counts of galaxies with 
15 < K < 25, and the inferred evolution o f the stellar 
mass function. In iPaper j we showed that the ! Bower et al.! 
model successfully reproduces the observed number counts 
of EROs, without having to adjust any of its parameters. 
The dust extinction is calculated self-consistently rather 
than being put in by hand (see Cole et al. 2000). We also 
predicted that EROs should be the most massive and bright- 
est of K-selected galaxies at 1 < z < 2, with a population 
dominated by old, passively evolving galaxies. Our predic- 
tions broadly agree with observations of EROs. We found 
that AGN feedback is a key ingredient for understanding 
the evolution of the most massive galaxies which were in 
place at 2; ~ 1. IFontanot fc Monaco! note that the modelling 
of stellar populations also plays a significant role in s haping 
the predictions for EROs (see also lTonini et al]|2009! ). 

In this paper we study the halo o ccupation distr i bution 
and clustering of EROs. Previously, I Quo fc Whit3 (|2009! ') 
reproduced the clustering of a similarly selected galaxy 



population, D RGs, using a different semi-analytical model. 
IGuo fc White! restricted their study to a single magnitude 
limited sample of galaxies. Here we will explore the vari- 
ation of ERO clustering with colour, magnitude, redshift, 
mode of star formation and also investigate how the cluster- 
ing changes from real to redshift space, trying to understand 
the physical processes behind each change. 

In this paper we select EROs using either their {R — K) 
or (i — K) colours. Colour cuts using (i — K) and also {I — K) 
appear to be more effective than (R—K) colours at removing 
"low" redshift gala xies with z < 0.8 and isolating galaxie s 
at higher redshif ts (jConselice et al.!l2008! : iKong et al.il2009l '). 
iKim et ahl (|201ll ) found that the clustering of EROs selected 
by their (R — K) or (i — K) colours is different, particularly 
at larger pair separations, where (i — K) selected galaxies 
have larger clustering amplitudes. 

The paper is organised as follows. In Section[2l we sum- 
marise the key features of iBower et al.l galaxy formation 
model. The predictions for the halo occupation distribution 
of EROs are presented in Section (3] The predictions for the 
spatial clustering are discussed in SectionUland those for the 
angular clustering in Section [5| Conclusions can be found in 
Section (H] 



2 GALAXY FORMATION MODEL 

We predict the clustering of EROs in a ACDM universe us- 
ing the GAL FORM semi-analyt ical galaxy formation model 
developed bv lCole et al.l l|2000l 'l. Semi-analytical models use 
simple, physically motivated recipes and rules to follow the 
fate of baryons in a universe in which structur e grows hier- 
archically through gravitational instability fsee lBaughll2006! . 
for an overview of hierarchical galaxy formation models). 

GALFORM follows the main processes which shape the 
formation and evolution of galaxies. These include: (i) the 
collapse and merging of dark matter haloes; (ii) the shock- 
heating and radiative cooling of gas inside dark matter 
haloes, leading to the formation of galaxy discs; (iii) qui- 
escent star formation in galaxy discs; (iv) feedback from su- 
pernovae, from active galactic nuclei (AGN) and from pho- 
toionization of the intergalactic medium (IGM); (v) chem- 
ical enrichment of the stars and gas; (vi) galaxy mergers 
driven by dynamical friction within common dark matter 
haloes, leading to the formation of stellar spheroids, which 
also may trigger bursts of star formation. The end product 
of the calculations is a prediction for the number and prop- 
erties of galaxies that reside within dark matter haloes of 

different masses. 

In this paper we focus our attention on the I Bower et al.l 
(j2006l ') variant of GALFORM. In lPaperl we showed that this 
mo del reproduce s the o bserved numbers of EROs, whereas 
the iBaugh et al.l (120051) model und erpredicted the counts 
3 iBower et al. model is implemented 

an N- 



of these galaxies. The iBower et al 



in the Millennium Simulation (jSpringel et al. 
body simulation with about 10^° particles, each with a mass 
of 8.6x lO*/i"^M0, in a box of side 500/i"^Mpc. The galaxy 
population is output at selected snapshots in the Millennium 
Simulation. The outputs are evenly spread in the logarithm 
of the expansion factor and so do not correspond to round 
numbers in redshift, as will become apparent when we plot 



© 0000 RAS, MNRAS 000, 000-000 



EROs clustering in a hierarchical universe 3 



our predictions (e.g. Fig. [T] shows predictions at 2 = 2.1 

rather than z = 2). 

Key features of the lBower et al.l model include (i) a time 
scale for quiescent star formation that scales with the dy- 
namical time of the disk and which therefo re changes signifi- 
cantly with redshift fsee lLagos et al.l|2010l . for a study of dif- 
ferent star formation laws in quiescent galaxies), (ii) bursts 
of star formation occur due to both galaxy mergers and when 
disks become dynamically unstab le, and (iii) feedbac k from 
both supernovae and AGN (see iBenson et al.l |2003| . for a 
discussion of the effect that feedback has on the luminosity 
function of galaxies). The onset of the AGN suppression of 
the cooling flow is governed by a comparison of the cooling 
time of the gas with the free-fall time for the gas to reach 
the centre of the halo. Cooling is suppressed in quasi-static 
hot haloes if the luminosity released by material accreted 
onto a central sup ermassive black hole balances the cool- 
ing luminosity fsee iFanidakis et al. I l201ll. for a full d escrip- 



tion of black hole growth in the model) . iBower et al.l adopt 
the cosmological par ameters of the Millennium Simulation 
l|Springel et al.ll2005l ). which are in broad agreement with 
constraints from measurements of the cosmic microwave 
background radiatio n and large scale galaxy clustering (e.g. 
ISanchez et a"l]|2009l ): fio = 0.25, Ao = 0.75, fit = 0.045, 
erg = 0.9 and h = 0.73, such that th e Hubble cons tant today 
is Ho — 100 h kms~^Mpc~^. The iBower et all model pa- 
rameters were fixed with reference to a subset of the avail- 
able observations of galaxies, mostly at low redshift (see 
iBower et aDl20ld . for a discussion of parameter fitting) . This 
model successfully reproduces the inferred stellar mass func- 
tion up to 2: = 4.5. For further det ails we refer the rea der 
both to IPaper j in this series and to lBower et~al] (|2006h . 

Here we extract p redictions for the clustering of EROs 
usmg theH ower et aU model, without adjusting any of its 
parameters. 

The bands used here correspond to the R band from 
SUBARU, the i band from the SDSS and the K band from 
UKIRT, with effective wavelengths of 0.65 fim, 0.75 fim and 
2.2 jj,m, respectively. All magnitudes used in this paper are 
in the Vega system, unless otherwise specified. The Vega-AB 
offsets are —0.196, —0.355 and —1.87 magnitudes for the R, 
i and K bands, respectively. 



3 THE HALO OCCUPATION DISTRIBUTION 
OF EROS AND THE ABUNDANCE OF 
THEIR HOST HALOES 

In this section we explore first the halo occupation distri- 
bution (HOD) of EROs and then their HOD-weighted halo 
mass function, which helps to relate the HOD to the clus- 
tering predictions. 



3.1 The HOD of EROs 

The spatial clustering of galaxies can be quantified by the 
correlation function, i.e. the excess or deficit of the num- 
ber of galaxy pairs at a given separation with respect to 
a random distribution. The correlation function of galaxies 
can be derived from the halo occupation distribution once 
the mass function and cl ustering of the haloes are known 
(e.g. iBenson et aU bOOObl ) . The steps relating the HOD to 
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Figure 1. The HOD at z = 2.1 of galaxies with K < 20 (solid 
black line) and of EROs with the same K limit and (R — K) > 
5, 6, 7 (dashed, dash-dotted and dash-triple dotted lines lines re- 
spectively, as indicated by the key) and (i — K) > 4.5 (dotted 
line). 



the effective bias of a galaxy sample have been reviewed by 
iKim et all (|2009[ ). 

The HOD gives, as a function of halo mass, the mean 
number of galaxies per halo, {N)m, which pass a partic- 
ular observational selection. Semi-analytical models natu- 
rally predict the form of the HOD, which is determined by 
the interplay between a range of physical processes which 
can change with redshift. The form of the HOD can be un- 
derstood in terms of the sum of the contribution of cen- 
galaxies, which typically is assumed to resemble a 
step function, and that of satellite galaxies, which approxi- 
mate ly follows a power law l|Berlind et al ] |2003l : IZheng et al.l 
l2005i ). This implies that, in general, {N)m has a cutoff 
at low masses and then a slow rise or pla teau followed 
by a po wer law with increasing ha lo mass l|Benson et al.l 
2000a||bt iBerhnd fc Weinberg! I2OO2I ). However, the HOD 
plat eau has been found to be less evident for brighter galax- 
ies (jZehavi et al.l l2010l ) and the same is expected to hap- 
pen for an older p opulation of galaxies (jBerlind et al.ll20C)3l : 
IZheng et al .120051 . note, however, that these theoretical stud- 
ies were developed before the inclusion of AGN feedback in 
galaxy formation models). 

In Fig. [1] the predicted HOD of EROs with K < 20 
at 2 = 2.1 is compared with that of all galaxies with the 
same magnitude limit, but without any colour cut. At this 
redshift, the HOD of K-selected galaxies and of EROs with 
either (R—K) > 5 or (i —K) > 4 .5 differ only in the low mass 
cutoff. As we found in iPaper^ this shows that EROs with 



^ In hierarchical models the central galaxy is defined as the most 
massive one within a dark matter halo. This galax y is, generally, 



situate d close to the centre of mass of the halo (|Benson et al.l 
l2000al lbl: iBerlind et aIll2003^ . 
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Figure 2. Solid lines show the HOD of all EROs with K <20 at 
z = 2.1. The dotted (dashed) lines show the predicted HOD for 
satellite (central) EROs. Different panels show the HOD predic- 
tions for EROs selected with different colour cuts, as labelled. 



K < 20 tend not to be hosted by haloes with low masses 
(Mhaio > 4 X lO"/i"^M0 in this case). At z = 2.1, the cut 
in apparent magnitude oi K — 20, corresponds to an abso- 
lute magnitude of Mk — 5\ogh ~ —25.5 (observer frame). 
According to the model, this value is close to L* at 2; = 2.1. 
Thus, Fig. [1] s hows the HOD of ~ L* and bri ghter galax- 
ies. In lPaper] it was shown that in the lBowei~et al. model, 
EROs account for most of the bright end of the K-band lu- 
minosity function at 1 < z < 2.5, explaining the similar 
HODs predicted for EROs and for all K-selected galaxies at 
2; = 2.1. At lower redshifts in the ERO range, moving toward 
z = 1, the difference between the ERO HOD and the overall 
galaxy HOD becomes larger. In this case, we are picking up 
less luminous galaxies and a smaller fraction of these satisfy 
the colour cut to be classified as EROs. 

Fig. [T] shows that the HOD of aU galaxies with K < 20 
at z = 2.1 is predicted to have a clear cutoff at low masses 
followed by an approximately power law rise. Below Mhaio ~ 
10^^h~^ Mq, the probability of hosting an ERO is close to 
zero. Nevertheless, as we shall see in the next subsection, 
such haloes where the HOD is rising from z ero can have a bi g 
influence on the predicted clustering (see lKim et al.ll2009l ). 
To model this transition region accurately it is necessary to 
generate many examples of merger histories of such haloes, 
to accurately quantify the fraction that host EROs. 

Fig. [1] shows that the HOD of all galaxies with K < 20 
at 2 = 2.1 flattens out at a value of {N}m below unity, 
in contradiction with the expectation for the generic form 
of the HOD. This implies that central galaxies with these 
characteristics are not expected to be found in every halo. 
The same trend is found for EROs. Fig. [2] shows that for 
EROs at z = 2.1, the HOD of central galaxies only reaches 
unity for haloes with M > 10^'^h~^ Mq, except in the case 
of EROs with {R — K) > 7, for which a mean occupation of 
unity is never reached. 

To understand this prediction, we examine the trend be- 



tween the K-band luminosity and the host halo mass. If there 
was a well deflned trend between the luminosity of a central 
galaxy and the mass of its host halo, then the central galax- 
ies residing in haloes above a given mass would typically be 
brighter than some limiting magnitude. We say "typically" 
to allow for scatter in the relation b etween galax y luminos- 
ity and halo mass. However, in the lBower et al.l model, we 
do not flnd such a correlation. For low halo masses, those 
in which AGN suppression of gas cooling is not effective, a 
monotonic relation between central galaxy luminosity and 
halo mass is predicted, with sufficiently small scatter to 
make the trend clear. This relation flattens and the scatter 
increases substantially beyond the mass in which gas cooling 
is first suppressed by accretion onto the central supermas- 
sive black hole. A central galaxy of a particular luminosity 
can be found in a broad range of halo masses, covering sev- 
eral decades in mass. For models including AGN feedback, it 
is no longer the case that the brightest central galaxies will 
be found exclusively in the most massive haloes. Instead, 
the tails of the luminosity-halo mass relation play a critical 
role. We find that for the brightest galaxies, only a fraction 
of haloes of any mass have formation histories which can 
produce such an object, hence {N)m ^ 1. 

Fig. [1] shows that the minimum mass. Mi , required for 
a halo to host two K-selected galaxies is a factor of ~ 11 
higher than that required to host only one, which is labelled 
Mmin- This is also the case for EROs except for those with 
{R — K) > 7. In this case, haloes are predicted to host at 
most one ERO with {R - K) > 7 and thus, {N)m < 2 for 
all masses. The similarity in the factor found for EROs and 
K-selected galaxies indicates that our result is not peculiar 
to the additional colour selection associated with an ERO 
sample, but appears to be a characteristic of bright galaxies. 
Indeed, we note that a similar ratio Mi/Mmin ~ 11 has 
been suggested for local galaxies br ighter than I . IL*, while 
Ml /Mmin ~ 17 for fainter galaxies (Zehavi et al.ll2010l ). 

Fig. [2] compares the predicted HOD for EROs selected 
with different colour cuts, separating the contributions from 
central and satellite galaxies. We note in Fig. [2]that the step 
in the HOD is less pronounced for redder central EROs. As 
an example, we find at 2 = 2.1 an average of 0.1 EROs 
with {R- K) > 5 and K < 20 per halo with 10^^/i"^Mq, 
while fewer than one in a thousand halo es host E ROs with 
{R — K) > 7 and K <.20 sX this mass. In lPaper^ , we found 
that redder cuts select brighter and older EROs. Thus, on 
the one hand, the variation of the HOD with colour pre- 
dicted here is in agreemen t with a (N )m shape dominated 
by br i ght and old galax i es (iBerlind et al . 2003; Zhe ng et al.l 
I2OO5I : ikim et all l2009l : IZehavi et al.l l2O10i ). On the other 
hand, Fig. [2] shows that the split between central and satel- 
lite galaxies varies with the colour cut used to define the 
ERO sample. The fraction of satellite galaxies increases for 
the redder cuts. The difference in the number of satellite 
EROs causes the shape of the halo occupation distribution 
to vary from a clear step function plus a power law for EROs 
with (R—K) > 5 or 6, to a distribution close to a pure power 
law for EROs with {R- K) > 7. In this case, the HOD for 
central galaxies does not get close to unity before it is dom- 
inated by the satellite HOD. 
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Figure 3. The mass function of all dark halos (solid lines) and the 
HOD-woighted halo mass function for different ERO selections. 
As in Fig.[T]dotted lines correspond to the mass function weighted 
by EROs with {i - K) > 4.5, dashed lines to {R- K)> 5, dash- 
dotted to (R- K) > % and dash-triple dotted to (R - K) > 7. 
Left-hand panels show the prediction for EROs with K < 19 and 
right-hand ones for K < 20. Top row corresponds to z = 1.1, the 
middle to 2 = 1.5 and the bottom row to z = 2.1, as labelled. 



3.2 The HOD-weighted halo mass function 

The HOD by itself does not contain any explicit informa- 
tion about the number density of dark matter haloes. For 
a given redshift and assuming an initial power spectrum of 
matter fluctuations, the number density of halos as a func- 
tion of mass, i. e, the halo mass function , depends only on 
the cosmology. iPress fc Schechtal l|l974 ) argued that for a 
Gaussian density field the abundance of haloes more massive 
than a characteristic mass is expected to decline exponen- 
tially, with a power law like behaviour at lower masses, a 
prediction confirme d by numerous N-body simulations (e.g. 
I Jenkins et aLllioOll ). Fig. |3] shows as solid black lines the 
mass function of dark halos at different redshifts predicted 
by the Millennium Simulation. 

We plot in Fig. |3] the halo mass function multiplied 
by the mean number of EROs per halo, where the EROs 
are defined with different colour cuts, {i — K) > 4.5 and 
{R — K) > 5, 6, 7, and in two magnitud e ranges, K < 19 
and K < 20. As anticipated in IPaper j we see from Fig. 



[3] that EROs can be found in halos with a wide range of 
masses. In the lowest mass bins we see a sharp cutoff driven 
by the shape of the HOD. 

We note that a colour cut not shown in Fig. |3] but that 
is widely used for observationally selecting EROs is {R — 
K) > 5.3. At z = 1.1, 1.5 and 2.1, the HOD for EROs with 
{i — K) > 4.5 is practically the same as that for EROs with 
(R-K) > 5.3. 

Fig. |3] shows how ai z = 1.5 and z — 2.1 the weighted 
mass functions for EROs with (i — K) > 4.5 and (R—K) > 5 
practically overlap, while they differ at z = 1.1. For the 
two magnitude limits shown in Fig. [3l the predicted redshift 
distributions of EROs with either (i — K) > 4.5 or [R—K) > 
5 have medians around z — 1.5. However, the low redshift 
tail cuts off more sharply for EROs with {i — K) > 4.5 and, 
therefore, these EROs are expected to contribute less to the 
mass function at z — 1.1, explaining the difference seen in 
Fig.El 

The mass of haloes such that (A'^) a/ = 1 for EROs is 
given by the mass at which the HOD-weighted halo mass 
function has the same amplitude as the halo mass function. 
In Fig. [3] this happens when a line crosses the dark mat- 
ter halo mass function (solid line). More massive halos will 
contain, on average, more than one ERO. In Fig. |3]we see 
that for EROs with {R - K) > 5 and K < 19 the mass at 
which their (TV) a/ = 1 increases by a factor of ~ 2.5 from 
z = 1.1 to z — 1.5. Fig. |3] shows again that EROs only 
inhibit halos more massive than a certain threshold value, 
which depends on redshift, apparent magnitude limit and on 
colour. Both the mass at which (N) m ~ 1 and the minimum 
mass for a halo to host an ERO increase with redshift. This 
trend is related to the fixed apparent magnitude limit which 
means that we are looking at intrinsically brighter galaxies 
at higher redshifts. 

In contrast, the median mass of haloes hosting EROs 
increases with reducing redshift. For EROs with (R—K) > 5 
and K < 19 the median mass of their host haloes increases 
from 1.8 X 10^^ h'^Me at z = 2.1, to 4.4 x 10^^ h'^Me at 
z = 1.1. This increase is close to that for the dark matter 
M*. Thus, this tendency is related to the growth of dark 
matter haloes with time and it is seen also for EROs with 
(i — K) > 4.5 and for K-selected galaxies. 

We have found that brighter EROs are hosted by more 
massive haloes. In particular, at z = 1.5, halos more massive 
than ~ 2 X 10^^ Mq will contain, on average, at least one 
ERO with K < 19, while one or more EROs with K < 20 
will be present in halos with masses ~ 5 x 10^^ Mq. Our 
pr edicted values agree with those from the empirical model 
of iMoustakas fc Somerviii^ (|2002l ). 

In general, as shown by Fig. |31 the redder an ERO, the 
more massive is the typical host halo. In particular, all halos 
are predicted to host on average fewer than one ERO redder 
than (R-K) = 7. For EROs with K < 19 at z = 1.5, 
the median mass of their host haloes increases from 2.5 x 
10^^ /i-^Mq for EROs with {R-K) > 5, to 6.3x10^^ ft-^M© 
for EROs with {R - K) > 7. 

We have explored the HOD-weighted mass function sep- 
arating the contribution from satellite and central EROs, 
and find that at z = 1.1 central EROs are rare, i.e. they 
make a marginal contribution to the ERO-weighted mass 
function for massive haloes. The average number of central 
EROs at z — 1.1 is well below one per halo, even for halos 
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as massive as 10"/i"^Mq. This is not the case at redshifts 
z = 1.5 and z = 2.1. At these redshifts, we find that massive 
haloes are expected to host an ERO as their central galaxy 
(except at z = 2.1 and K < 19, in which case EROs are just 
rare). Thus, the increase from z = 1.5 to z = 1.1 in the halo 
mass where {N)m = 1 for EROs with (i — K) > 4.5, appears 
to be related to the fact that at lower redshifts the fraction 
of satellite EROs increases. 



4 THE TWO-POINT CORRELATION 
FUNCTION OF EROS 

In this section we present the predictio ns for the two -point 
correlation function of EROs, ^(r). The lBower et al] model 
is implemented in the Millennium simulation and thus, it 
can be used to directly predict the spatial distribution of 
galaxies and hence the two-point correlation function. In 
real-space, the Cartesian comoving co-ordinate of EROs 
within the simulation box are used to compute pair separa- 
tions. We calculate their two-point correlation function, ^, 
using a simple estimator (e.g. [Peebles 1980.): 



1 + ^ 



2DD 
nWdV 



(1) 



where DD is the number of distinct galaxy pairs with sepa- 
ration between r and r + dr measured from the simulation, 
allowing for periodic boundary conditions, n is the mean 
number density of galaxies, V is the total simulated volume, 
and dV^ is the volume of the spherical shell of radius r and 
thickness dr. The denominator of Eq. [T] corresponds to the 
average number of neighbours found in dV^ in the case of a 
Poisson distribution. 



4.1 The shape of the two-point correlation 
function 

The predicted real-space two-point correlation function at 
2; = 2.1 for galaxies with K < 20.2, including EROs, is plot- 
ted in Fig. [4] This magnitude limit has been chosen, to match 
that used in deep ob servational surveys jRoche et al.ll2002l : 
iMivazaki et al.|[2003l ). The predicted shape of the two-point 
correlation function for all K-selected galaxies and EROs is 
roughly consistent with a power law for pair separations in 
the range 0.03 < r < SO/i'^Mpc. 

Fig. [5] plots the ratio between the correlation function 
of the EROs and a power law, ~ (r/ro)^, of EROs with 
{R - K) > 5 and K < 20.2 at z = 2.1. This emphasizes 
any deviation from a power law of the real-space correlation 
function. The parameters of the power law ^flt are those that 
give the best fit to the real space correlation function of the 
EROs plotted in Fig.[5l ro = 7.0 h'^Mpc and 7 = -2.1. Fig. 
[S]shows that there is a change in slope at r ~ 1/i^^Mpc. We 
have quantified this change by fixing ro and then fitting ^ 
separately for r < l/i~^Mpc and for larger separations. We 
have found a change in slope IA7I < 0.5 (less than 25%) 
between large and small scales for the case of EROs selected 
with different colour and K-magnitude cuts and at redshifts 
z — 1, 1.5 and 2. 



^ Unless otherwise noted, hereafter we will refer to comoving 
co-ordinates simply as coordinates. 



O 




i-K)>4.5 

Open: 1 —halo term 
Filled: 2 — halo term 




Figure 4. The predicted real-space two-point correlation func- 
tion, ^, at 2: = 2.1, as a function of comoving separation. The 
upper panel shows the predicted ^ for all K-selected galaxies with 
K < 20.2 (dashed line), and for those which are also EROs, with 
(i — K) > 4.5 (triangles), (i? — K) > 5 (circles), 6 (squares), 7 
(pentagons). The bottom panel shows the predicted ^ for EROs 
with (i — K) > 4.5, distinguishing between the contributions from 
the one-halo term (open triangles) and the two-halo term (filled 
triangles; see text for details). In both panels we also show the 
predicted ^ of the dark matter (solid line). The dotted line shows 
^ = 1 and the shaded areas show the Icr Poisson errors derived 
using the number of pairs predicted by the model at a given sep- 
aration. 



Gravitational instability theory predicts the existence 
of an inflection point in the two-point correlation function 
of dark matter when a smooth initial power spectrum is as- 
sumed, such as a power law. This occurs at the transition 
from the linear to the non-linear regime, whic h is rela.ted to 
the one-halo to two-halo transi ti on (see e.g. Peebles 1980l : 
iGaztafiaga fc Juszkiewic3 I2001I : IScoccimarro et al.l 200in . 
From both semi-analytical models and smoothed particle 
hydrodynamics simulations, a similar change in slope is ex- 
pected to occur for the two-point correlation function of 
galaxies, for separations close to the dark matter correla- 
tion length, ro, i.e. the separation for which ^DM{ro) = 1. 
The strength of the inflection depends on the galaxy se- 
lection which determines the number of galaxies per halo 
and the balance between the one and two halo term. This 
transit i on is supported by observations at low redshift (e.g. 
Baughl 1 19961 : iGaztanaga fc Juszkiewic3 I2001I : IZehavi et al.l 



200 j ). We can see such behaviour in Figs. |4]and[5l Around 



l/i~^Mpc, the two-point correlation function measured for 
all K-selected galaxies and EROs has a small change of slope 
that makes the correlation function steeper on small scales. 
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Figure 5. The predicted ratio between the spatial correlation 
function of EROs, with {R - K) > h and K < 20.2 at z = 2.1, 
and a power law fit, ^fjt = {r/ra)'^ , where ro is 7.0/i~^Mpc, and 
7 is —2.1. The solid black (grey) line corresponds to the ratio 
using the real (redshift) space correlation functions. The dashed 
line is the redshift-space clustering for r > lh~^Mpc obtained by 
multiplying the clustering in real-space by the Kaiser factor (Eq. 
[21 i|4.1.4| l. The shaded areas show the Icr Poisson errors derived 
from the number of pairs at a given separation. 



On small scales, r < lh~ Mpc, ^ mainly measures the 
number of pairs o f galaxies within the same halo, i.e. the 
one-halo term (e.g. iBenson et al]|2000bl ). The one- halo term 
is dominated by the distribution of satellites within single 
haloes. Thus, a sample with a higher proportion of satel- 
lites will have a boost in the clustering on small scales. The 
distribution of satellite galaxies depends strongly on the se- 
lection criteria (magnitude limit, colour cut, redshift range, 
etc.). The lower panel of Fig. [4] shows separately the con- 
tribution of the one-halo and two-halo terms to the pre- 
dicted global clustering of EROs with {i — K) > 4.5. We 
note that the two predicted contributions for EROs with 
{R — K) > 5 are almost indistinguishable from those for 
EROs with (j - K) > 4.5. 

On larger scales, r > 1 ^~^Mpc, we are generally prob- 
ing the positions of galaxies in different halos, the two-halo 
term. We can see in the lower panel of Fig. U that the 
two-halo term of EROs has the same shape as the corre- 
lation function of the dark matter. On these scales, the two 
point correlation function of galaxies traces that of their 
host dark matter haloes scaled by their bias and the dark 
matter clustering is close to the prediction from linear the- 
ory on these scales. We measure the clustering of EROs 
with {R - K) > 5 a.nd K < 20.2 a,t z = 2.1, ^gg{r), to 
be boosted with respect to that of dark matter, ^dm, with 
a bias b ~ \/£,gg{r)/(,DM ~ 3 (at r = 8/i~^Mpc). This is 
consistent with these galaxies being hosted by haloes with 
Mhaio ~ 4 X lQ^^h~^ Mq, in agreement with the predictions 
shown in Fig. 13 

In the lower panel of Fig. [4]it is noticeable that the two- 
halo term for EROs does not go to zero at scales where the 
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Figure 6. Variation with redshift of the predicted real-space 
two-point correlation function for EROs with (H — K) > 5 and 
K < 20.2. Correlation functions at 2 = 1.1, 1.5, 2.1, 2.4 (dashed, 
dotted, solid and dash-dotted line, respectively) are shown di- 
vided by the correlation function predicted at z = 2.1. 



one-halo term dominates. For EROs there is a large satellite 
component. It is then possible for two haloes that are almost 
touching, to have satellite pairs at a smaller separation than 
that corresponding to the centres of the two haloes. 

4-1.1 Variation of clustering with colour 

Fig. [4] shows that at z = 2.1 galaxies in the ERO samples 
defined by the bluest cuts, with {R — K) > 5 or (i — K) > 
4.5, cluster in a very similar way. In line with the results 
described in [JS] the clustering of these bluest EROs is also 
seen to be very close to that of galaxies selected solely on 
the basis of their K-band apparent magnitude without a cut 
in colour. EROs with {R — K) > b or {i — K) > 4.5 account 
for most of the bright end of the luminosity function of K- 
selected galaxies at z = 2.1, which explains their similar 
clustering. 

Fig. [4] also shows that the reddest EROs display the 
strongest change in power law slope with pair separation. 
As discussed above, a larger proportion of satellite galaxies 
is found among the redder EROs, explaining the boost seen 
at small scales. 

4- 1-2 Variation of clustering with redshift 

Table [T] presents the real-space comoving correlation length, 
ro, of EROs with [R — K) > 5 for different redshifts and 
magnitude ranges. The correlation length is obtained here 
as the galaxy separation that gives ^(ro) = 1. 

In Table [T] we can see that ro does not change monoton- 
ically with redshift. This lack of a cle ar evolutionary tren d 
has also been found observationally bv I Brown et al.l (|2005h . 
A simil ar prediction was m ade from the theoretical predic- 
tions in lBaugh et al.l (|l999l ) for a sample with a fixed appar- 
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Table 1. The predicted real-space comoving correlation length, 
ro, for EROs with (i? — K) > 5 and different magnitude cuts, at 
2 = 1.1, 1.5, 2.1. Not enough EROs with K < 18.4 are found at 
2 = 2.1 to have an estimation of their clustering. 
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ent magnitude limit, so this trend is not necessarily peculiar 
to the colour selection. 

The variation of the clustering of EROs with redshift is 
small. Therefore, we have emphasized this variation by pre- 
senting in Fig. |6] the clustering of EROs with {R — K) > 5 
and K < 20.2 at different redshifts, z = 1.1, 1.5, 2.1 and 
2.4, divided by that ai z = 2.1. It is clear from Fig. [S] 
that the two-point correlation function of EROs in real-space 
presents a slightly different shape at z = 1.1 than at higher 
redshifts, showing a boost on small scales r < 0.5/i~^Mpc. 
This is related to the higher fraction of satellite galaxies se- 
lected at z = 1.1 as EROs, in comparison with the higher 
redshifts plotted, as reported in 321 

The evolution of the correlation function with redshift 
is important for interpreting the angular clustering of galax- 
ies. This evolution is often parametrised as a global, scale 
independent factor in redshift. Fig. [6] shows that the varia- 
tion of the clustering with redshift is different on large and 
small scales. This makes it difficult to assign a unique red- 
shift power law, 1/(1 + 2)", to describe this evolution, so 
such "deprojections" of angular clustering are at best ap- 
proximate. 
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Figure 7. The predicted spatial two-point correlation function of 
EROs, with (R- K) > 5, brighter than K < 18.4 (circles), 19.2 
(squares), 20.2 (triangles) at 2 = 1.1. The real-space prediction is 
shown in the upper panel and the redshift-space in the bottom. 
In both cases, the shaded areas show the Icr Poisson errors, de- 
rived using the number of pairs predicted at a given comoving 
separation. 



4-1.3 Variation of clustering with apparent magnitude 

Fig. [7] shows the predicted two-point correlation function 
at z = 1.1 in both real and redshift-space for EROs with 
{R~ K)>^ &nA K < 18.4, 19.2 and 20.2. At this redshift, 
the model predicts the existence of EROs in sufficient num- 
bers for each of these magnitude cuts to allow as to compare 
the clustering characteristics. We find that the clustering of 
EROs depends weakly on luminosity. Fainter EROs are pre- 
dicted to be slightly more clustered than brighter ones, due 
to a higher fraction of satellite galaxies in massive haloes. 
Fig. [7] shows that the strongest variation occurs on small 
scales. Table [1] also shows a weak trend for fainter galaxies 
to have a higher comoving correlation length, ro. 

Table [T] lists the predicted ro for EROs selected in 
magnitude ranges similar to various observational stud- 
ies that inferred the correlation length by fi tting models 
to the measured angular correlation functions. lOaddi et al.l 
l|200lh analysed 400 EROs with Ks < 19.2, [R - Ks) > 5 
and a median redshift z ~ 1.2 in a 701 arcmin^ field and 
estimated that their correlation length is not less than 
ro — 8/t~^Mp c, with a most pr obable estimate of ro ~ 
12±3fe~^Mpc. lRoche et all (|2002l 'l inferred that EROs with 
K < 20.2, (7? - Tf) > 5 and 1 < Zphoto < 3 (158 EROs 
in a 81.5arcmin^ field), h ave a correlat i on len gth between 
10/i~^Mpc and ISfc'^Mpc. ISrown et all (120051 ) inferred the 



correlation length of EROs with K < 18.4, {R~ K) > 5 and 
a median photometric redshift of 2 ~ 1.2 (671 EROs in a 
3529arcmin^ field), to be 9.7 ± 1.0 /i^^Mpc. 

Table [1] shows that the predicted correlation length for 
EROs is consistently lower than that estimated from the 
observations. However, the predicted correlation length for 
EROs with K < 19.2 at z = 1 .1 is above the lower limit 
estimated bv lDaddi et"al] (|200ll ) and within the 2a range of 
their best estimate of ro. The comparison with correlation 
lengths deduced from observations is not ideal since a num- 
ber of assumptions have to be made to infer this quantity 
from the actual observables. In particular, all the observa- 
tional studies mentioned above assumed a power law for the 
spatial two-point correlation function, ^, with a fixed index 
7 « — 1.8. The predicted ^ is best fitted by a power law with 
an index 7 « —2.1 at z = 1.1 and at 2 = 2.1, as shown at the 
beginning of H4.ll Thus, if this predicted value was adopted 
for the observational estimation, lower correlation lengths 
would be obtained, resulting in a better agreement between 
the model and observations. A direct comparison with the 
observed two-point angular correlation function will be car- 
ried out in [J5l 
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4-l--i Clustering in real and redshift- space 

To mimic a spectroscopic survey in which radial positions 
of galaxies are inferred from their redshifts, galaxy positions 
are perturbed along one of the axes by their peculiar ve- 
locities, scaled by the appropriate value of the Hubble pa- 
rameter. The impact of these peculiar motions depends on 
the scale. Figs. [5] and [7| show that while the slight change of 
slope around l/i~^Mpc in the correlation function of EROs 
is clear in real space, this is smeared out when including 
redshift-space distortions. 

Fig. [5] emphasizes the differences between the correla- 
tion functions in real and redshift-space by plotting the cor- 
relation functions divided by the power law that best fits 
the real-space clustering of EROs. Both Figs. [5] and [7] show 
that the difference between real and redshift-space cluster- 
ing depends on scale. 

On small scales, r < l/i~^Mpc, the clustering in 
redshift-space is significantly lower than that in real space. 
On these scales we are generally considering ga laxies within 
the same halo. As can be seen in Fig. [3l the I Bower et al.l 
model predicts the existence of more than one ERO in halos 
of high enough mass, except for the most extreme EROs with 
(f? — K) > 7. The peculiar motions of EROs within a halo 
cause an apparent stretching of the structure in redshift- 
space, diluting the number of ERO pairs at small sepa- 
rations. Thus, the correlation function signal measured in 
redshift-space is smaller than that in real-space on these 
scales. 

On larger scales, 1 < r < 30/i^^Mpc, bulk motions of 
galaxies cause the galaxy distribution to appear squashed 
along the line of sight in redshift-space, enhancing the clus- 
tering amplitude. Fig. [S] shows that, on this range of scales, 
the boost measured for the redshift-space correlation func- 
tion is in reasonable agreement with th e value expected ac- 
cording to the formalism developed by ' Kaiseil (|l987l V This 
relates the correlation function in redshift-space, £,gg{s), to 
that in real-space, S,gg{r), through a factor, /, ^gg{s) — 
f£,gg{r). This factor is a function of the bias, b, and the 
matter content of the Universe, ilmi 



/ = 1 + 



3 b 



+ 



5\ b J 



(2) 



In fact, the redshift-space clustering is only expected to 
be accurately described by the "Kaiser factor" on very 
large scales, since li near perturbation theory is assumed 
l| Jennings et aLlbOllI ). Tr aditiona,lly 7' = 0.6 but 7' = 0.55 
is a better approximation (|Linderll2005l ) and is the value we 
have used here. 



4.2 Clustering of quiescent and star forming 
EROs 

Observationally, EROs exhibit a mixture of spectral prop- 
erties which are related to their star formation histories and 
dust content. In this section we analyse the differences in the 
predicted clustering between EROs that are actively forming 
stars and those that are passively evolving. 

Fig. [8] shows the clustering of EROs with {R- K) > 5 
and K < 20.2 at z = 1.5, divided into burst (galaxies with 
an ongoing burst of star formation), post-burst (last burst 
happened within the past 1 Gyr) and quiescent galaxies (all 
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Figure 8. The predicted real-space two-point correlation function 
for EROs with {R - K) > 5 and K < 20.2 at z = 1.5, separated 
into burst (squares), post-burst (triangles) and quiescent (circles) 
galaxies (see text in i|4.2l for details). 



the others). As described in iPaper ll . this classification is 
made by considering the lookback time to the start of the 
last burst of star formation experienced by the galaxy. In the 
I Bower et al.l model bursts are triggered by galaxy mergers 
and by disks becoming dynamically unstable. 

The higher clustering amplitude predicted in real-space 
for burst EROs, ro = 11.7/i~^Mpc, is noticeable compared 
with that of quiescent ones, for which ro = 7.3/i~^Mpc, and 
for post-burst ones, with ro = 6.7/i~^Mpc. However, Poisson 
errors dominate the burst ERO clustering estimate, due to 
the much lower space density of these EROs compared with 
the other two samples. We can see in Fig. |8]how the correla- 
tion functions for quiescent and post-burst EROs are similar 
for separations larger than r ~ l/i~^Mpc, with slopes differ- 
ing by IA7I < 0.2 on small scales. The same tendencies are 
found i n redshif t-space. 

In iPaper j we explored the predicted mix of quies- 
cent and active EROs using different classifications. In par- 
ticular, we used the colour-colour diagram proposed by 
iPozzetti fc Mannuccil (|2000t ). the specific star formation rate 
of EROs, adopting SFR/M^ = 10""yr"^ as the bound- 
arjjfl, and the classification described above, but consider- 
ing as active both post-burst and burst EROs. We found 
that the number counts of quiescent EROs agreed among 
the different classifications. 75% of the quiescent EROs 
in Fig. mhave SFR/K'h < lO'^yr'^ We find that the 
EROs with SFR/Ph < 10"^Vr"^ have a correlation length 
ro = 7.7/i~^Mpc and that, globally, their two-point correla- 
tion function is very close to that for quiescent EROs. The 



-1 

however there are very few EROs with SFR/Mt > 10 ~^°yr~ ^ 



At 2 = 1.5, the inverse of the age of the Universe is 10 yr 



The specific star formation boundary value was chosen in lPape^^ 
taking into account both the predicted distribution and number 
counts of EROs. 
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correlation lenght of EROs with SFR/M, > 10~"yr~^ is 
predicted to be ro — 5.7/i^^Mpc, slightly lower than that 
for post-burst EROs. 

The active EROs, post-burst plus burst EROs, are domi- 
nated in our model by the post-burst galaxies. The predicted 
clustering of active EROs practically overlaps with that for 
post-burst EROs only, having the same predicted correlation 
length ro = 6.7/i"^Mpc. 

A few observational studies have tried to obtain the 
real-space correlation length differentiating between quies- 
cent EROs and du sty star forming or active EROs. Among 
these are those by iDaddi et al.l l|2002l ) and iMivazaki et al.l 
l|2003l) studies, which inferred the correlation length of EROs 
from the measured angular clustering. 

Miyazaki et al. (200l) observed EROs with {R~Ks)> 
5 and Ks < 20.2, dividing the sample into quiescent and 
star forming galaxies by fitting observed broad band colours 
to synthetic spectral energy distributions. Without taking 
into account the photometric redshift errors, IMivazaki et al.l 
found the correlation functions of quiescent and active EROs 
to be ro = 11 ± 1 h'^Mpc and ro = 12 ± 2 fc'^Mpc, respec- 
tively. In agreement with the observations of lMivazaki et aU 
we find that quiescent and active EROs have similar corre- 
lation lengths. Taking into account the photometric redshift 
errors, IMivazaki et al.l revised their estimate of the correla- 
tion length for quiescent EROs to be ro = 8ft" ^Mpc. This 
value is comparable to but higher than the predicted cor- 
relation length of quiescent EROs with {R — K) > 5 and 
K < 20.2 at z = 1.5. 

iDaddi et all (|2002l ') observed EROs with {R~ R's) >5 
and Ks < 19.2, s plitting the sample using spectral fea- 
tures. iDaddi et al.l estimated the correlation length to be 



ro 



< 2.5 h ^Mpc for star forming EROs and 5.5 < 



ro 



16/i~^Mpc for quiescent EROs. Our predicted ro values 
for quiescen t EROs agree within errors with those from 
iDaddi et alj. The correlati o n fun ction for quiescent EROs 
inf erred by Mivazaki et al.l (j2003l ') agrees with the estimate 
bv IPaddi et a l. (2002), though they used different magni- 
tude cuts. However, these two observational studies esti- 
mated very different correlation lengths for star forming 
EROs. Sample variance could be responsible, at least partly, 
for the diffe rence between the observational studies, since 
IPaddi et al.l covered a field 7 times larger than that used 



the IMivazaki et al.l study. Another issue that will increase 
the disc repancy is the differe nt depths of the two obser- 
vations. ICimatti et aD (|2002ch observed that star forming 
EROs appear at slightly higher redshifts than old EROs. 
This implies that surveys sampling only bright galaxies will 
contain a lower percentage of star forming EROs. The need 
for more fields to be explored observationally is clear, in or- 
der to have a better knowledge of the origin of such different 
ro values for active EROs. 



5 THE TWO-POINT ANGULAR 

CORRELATION FUNCTION OF EROS 

In order to make a direct comparison with current measure- 
ments of clustering, we study here the predicted angular 
correlation function ijj{6) for EROs, obtained from our pre- 
vious estimate of the two-point spatial correlation function 
in real space and the redshift distribution of model EROs. 



The Limber equation l|Limbeilll95'i ) relates the spatial 
correlation function ^ to the angular correlation function cj, 
under the assumptions; 

(i) We can apply a sample selection and measure the cor- 
relation function for that full sample, without any concern 
for dependencies of clustering on i ntrinsic galaxy properties 
which vary within the sample (e.g. |Peebleslll980l ). 

(ii) The selection function does not vary over the pair 
separations on which we can measure a signal for ^ (see e.g. 
ISimon|[2007l ). For most surveys this is a reasonable assump- 
tion, since the signal for ^ will be measurable only for those 
pairs of galaxies such that their pair separations |ri — r2| are 
much smaller than their comoving distances, ri, r2. 

The above approximations are both reasonable for EROs. 
Therefore, for a fiat Universe we can calculate the angular 
correlation function, oj, for pairs of EROs separated by an 
angle 6 by: 



2/o°°[^]'if (/o"d<(r,2,.))d. 



(3) 



where dN/dz is the redshift distribution of the surveyed 
galaxies. For a pair of galaxies at comoving distances ri and 
r2: It = ri — r2, r = (ri -|-r2)/2 and the comoving separation 
between them is approximated by ri2 = \/u^ -\- r^vj^, with 
zu'^ = 2(cosS — 1). The variation of redshift with comoving 
distance dz/dr = (_ffo/c)-\/f2o(l + z)^ -I- (1 — f2o), where Ho 
is the Hubble constant now and c is the speed of light. When 
integrating the two-point correlation function at very large 
scales, beyond those modelled accurately within the simu- 
lation volume, we assume that on these scales ^ is given 
by a scaled dark matter two-point correlation function cal- 
culated from the linear initial power spectrum of density 
fiuc tuations, where the scaling is the linear bias factor, b^ 
(see lOrsi et"al]|2008l ). 



5.1 EROs selected on (R-K) colour 

Fig.djshows the predicted w{6) for EROs with (R-K) > 5 
and different apparent magnitude cuts, compared with ob- 
servations. The t riangles in F i g. sh ow the clustering of 
EROs studied by iDaddi et all j2000bl ') in a field, Daddi-F, 
which covers 701 arcmin'^ and is 85% complete to Ks < 18.8. 
The upside down triangles in Fig. [§ ] show the clustering of 
EROs studied bv lKong et al.l (|2006i ) in both a subsection of 
the Daddi-F field of 600 arcmin^ and in the Deep3a-F field 
with an area of 320arcmin'^. The s quares in Fig . sho w the 
clustering of EROs measured by iRoche et al.l (|2002| ) in a 
81.5arcmin^ field. The pentagons in Fig.[9] show the c l uster - 
ing of EROs with K < 17.9 measured bv iBrown et all (|2005t ) 
in a 0.98 deg^ field; their error estimate uses a Gaussian ap- 
proximation to the covariance matrix and, thu s, includes 
the eff ects of sample variance and shot noise. iKim et al] 
l|201ll ) measured the clustering of EROs with {r — K) > 5 
and K < 18.8 observed in the SA22 field with an area of 



For open cosmologies s ee the general expression given by 



tor open cosmologies S ( 
iBaugh fc Efstathioul ^99^ 
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Figure 9. The two-point angular correlation function of EROs 
with {R - K) > b and brighter than K = 18, 18.8, 20, from 
top to bottom. The black solid lines show the predicted angu- 
lar correlation. The top panel also contains a dotted line with 
the clustering predicted for EROs with K < 17.9. The mid- 
dle panel includes a dashed line showing the predicted angu- 
lar clustering for EROs with (r - K) > 5 and K < 18.8. 
In each panel t here are obs e rvation al results fo r the corre- 
spond ing EROs: iDaddi et all ll2000bl) (tria ngles), iRoche et all 



(I2OO2I') (filled squ ares), iBrown et al.r ||2005|) (filled pentag ons) , 
iKong et ai] |20oi) (upside down triangles) . iKim et al.lll201ll ) (cir- 
cles) .ErrorbarsarePoisson error bars, except for the data points 
from lBrown et al.l | |2005|) , which include the effects of sample vari- 
ance and shot noise. 



2.45 deg^, one of the lar gest fields used t o date to mea- 
sure ERO clustering. The lKim et al.l (l201ll ) data points are 
shown as open circles in Fig. 

Both iDaddi et al.] (|20()0bl ) and iKong et al.1 (I2006D ob- 
served roughly the same field using different R-filters. Their 
estimat es of the clu s tering of EROs agree within the errors, 
though [Kong et al.l (|2006|) measure d a slightly higher clus- 
tering amplitude than Daddi et al.l 

The top panel in Fig. [9] shows the predicted angular 
clustering of EROs for both K < 18, solid line, and K < 
17.9, dotted line, to match the observational samples. The 
predicted difference is minimal, as expected for such a small 



change in magnitude. The top panel in Fig.[9]s hows that the 
model reproduces the clustering estimated bv lBrown et al.l 
({20051 ) for EROs wit h K < 17.9 ii i a fiel d 5 times larger 
than that observed bv lDaddi et al.l l|2000bl ) and lKong et al.l 
(2006) (Daddi-F field). 

The middle panel in Fig. shows the predicted angu- 
lar clustering of EROs with K < 1 8.8 compared with the 
estimates f r omlDaddi et al.1 (l2000bh , iKong et al.1 (|2006l ) and 
iKim et al] (|201ll ). As we found for brighter EROs, our pre- 
dicte d angular clu steri ng is lower than that estimated by 
both lOaddi et al.l and iKong et al.l (Deep3a-F field) . How- 
ever, it is clear from Fig. [5]that the model predictions match 
the clustering for EROs with K < 18.8 estimated by the 
UKIDSS team llKim et al.ll201ll) over an area 12 times larger 
than the Daddi-F field. The middle panel in Fig. [5] also 
shows the predicted angular clustering for EROs selecte d 
with (r — K) > 5, the same filters used bv lKim et all l|201lh . 
Both r and R filters are very close, and the difference in the 
predicted angular clustering is negligible. 

For EROs with {R - K) > 5 and K < 20, the bot- 
tom panel of Fig. |5] shows that the predicted clustering re- 
produ ces t he observational es timates from both lRoche et all 
(I2OO2I ) and lKong et~aLl (|2006l ) (Deep3a-F field). This is also 
true for EROs with K < 19.5. 

iKong et al.1 (I2OO6I ) found an increase in the amplitude 
of the angular correlation function with the brightness of 
EROs. The observed variation of clustering with magnitude 
is stronger than that predicted, even after taking into ac- 
count some small discrepancies due to the use of slightly 
different filters. This is not peculiar to EROs , since a simi- 
lar conclusion was reached bv lKim et al.l l|2010t ) , who studied 
bright galaxies at z = 0, finding that the clustering pre- 
dicted by semi-analytical mod els does no t vary with lumi- 
nosity as strongly as observed. iKim et al.1 l|2010l ) concluded 
that a stronger dependence of clustering strength on lumi- 
nosity could be obtained by including additional physical 
processes which infiuence the mass of satellite galaxies. This 
could have an implication for EROs as an appreciable frac- 
tion of our EROs are satellites. 

Nevertheless, once the sampling variance is taken into 
account, the model predictions do match the observed an- 
gular clustering reasonably well for 0.0006 < 6'(deg) < 0.6. 

5.2 EROs selected on (i-K) colour 

Accurate measurements of the clustering of EROs can only 
be achieved by d eep imaging of l arge areas. Currently, the 
UKIDSS survey l|Kim et al.ll201ll ) is one of the few surveys 
that has been able to measure the clustering of EROs beyond 
separations of 0.1 deg (which corresponds to ~ 2/i~^Mpc 
at 2 = 1.5), which probes the two-halo component of the 
correlation function. In this section we show the model pre- 
di ctions for EROs se lected in a similar way to that adopted 
bv lKim et al.l (I2OIII ). i.e. using the (i — K) > 4.5 colour cri- 
terion. This criterion appears to select galaxies with z > 1, 
with fewer lowe r redshift contam inants than the {R — K) 
colour criterion l|Kong et al. I I2OO9D . 



5.2.1 Number counts and redshift distribution 

The success of a model in reproducing the observed angular 
clustering of a certain type of galaxy is more robust if the 
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Figure 10. DifFerential number counts of EROs with (i — K) > 
4.5. The predicted number counts for all EROs are shown by a 
solid line. The prediction for "old" ( "dusty" ) EROs are presented 
with dotted (dashed) lines (see i|5.2.3l for details on t his di vision). 
The triangles show the observations from lKim et al.l 1I2OIII ) , which 
are complete up to < 18.8. The inset shows the redshift distri- 
bution for EROs with (j — K) > 4.5 and K < 18.8: predictions 
are in black and observations in grey. The areas of the histograms 
are normalised to unity. The median redshifts are shown by the 
vertical arrows. 



model also reproduces their number density. Here we explore 
the predicted number counts and redshift distributions for 
EROs selected by their (i — K) colours. We already showed 
the predicted abund ance and redshift distributions for (R — 
K) selected EROs in lPaper j . 

Fig. [To] shows the predicted differential number counts 
for E ROs with (i—K) > 4.5, compared with the observations 
from lKim et"aLl (|201ll ). These observations are complete for 
galaxies with K < 18.8 (see iKim et allboill . for further 
details). The model underpredicts the number counts at the 
bright end, but gives a good match from = 18 to the 
completeness limit. 

The inset in Fig. [10] shows the predicted redshift dis- 
tribution of EROs with {i - K) > 4.5 and K < 18.8. 
For comparison we have i ncluded, in grey, the redshift 
distribution estimated by iKim et al] (|201ll ). This esti- 
mate is based on photometric redshifts mea sured by the 
NEWFIRM Mediu m Band Survey ( Bramm er et al.l I2OO9I : 
Ivan Dokkum et al.l(2oTol '). The area sampled by the NEW- 
FIRM survey is rather small, 0.25 deg^, and so is subject to 
significant sample variance, with the result that distinctive 
features appear in the redshift distribution. The model pre- 
dicts a redshift distribution close to the observed one, but 
with a slightly higher median redshift z = 1.44, as opposed 
to a median of 2 = 1.31 for the observational data. 

The model predicts a population of EROs defined with 
{i — K) colours that is close enough in number and redshift 
distribution to the observations to make it worthwhile to 
continue calculating their clustering properties. For this pur- 
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Figure 11. The predicted two-point angular correlation func- 
tion of EROs with {i — K) > 4.5 and K < 18.3 in the top panel, 
K < 18.8 in the bottom panel. The solid black lines show the pre- 
diction from the model. The dashed black lines show the model 
predictio ns using the observationally estimated redshift distribu- 
tion from iKim et al. I 1I2OIII ). The grey lines in the bottom panel 
present the predicted contribution of the one halo (dotted line) 
and two ha l o term s (dash-dotted line). Observational results from 
iKim et al.l boill ) for the UKIDSS DXS SA22 field are shown as 
open circles. The filled circles pr esent the observational estimatio n 
for the UKIDSS Elais-Nl field l lKim et al. in preparationll201ll) . 
Observational data points are shown with Jackknife error bars. 



pose, we will follow the method described above for EROs 
selected based on their {R — K) colour. 



5.2.2 Angular clustering 

Fig. [TT] shows the predicted two-point angular correla- 
tion function of EROs with (i — K) > 4.5 compared 
with the observed angular correla tion function me asured 
within two UKIDSS fields: SA22 l|Kim et ahlboill) . open 
circles, with an area of 2.45 deg- ^ and Elais-Nl, filled circles 
l|Kim et al. in preparation! |2011^ with an area of 3.88 deg'^ 
(both quoted areas are after masking). The SA22 field was 
observed with both the CTIO and UKIRT telescopes, while 
the Elais-Nl field was observed with the Subaru and UKIRT 
telescopes. The observed correlation function in both fields is 
corrected for the integral constraint. For the Elais-Nl field 
the integral constraint is estimated to be 0.0073, which is 
comparable to the clustering signal at scales of « 0.7 deg. 
For the clustering i n both fields Jackknife e rrors have been 
estimated following ISawangwit et al] (|2009l ), in an attempt 
to take into account the sampling variance (|Norberg et al.l 
|2009i). 

We predict the angular clustering using Limber's equa- 
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tion (Eq. [3]), assuming either the predicted (solid hnes) 
or the observed (dashed hnes) redshift distribution for the 
galaxies under study. We can see from Fig. [TT]that the dif- 
ference between using the predicted and observed redshift 
distributions is negligible. The evolution with redshift of the 
spatial two-point correlation function for EROs is modest, 
which explains the small variation seen for the angular clus- 
tering when using a different r edshift distribution. 

From the lKim et al.l (|201ll ) observations, the correlation 
function of EROs selected with {i — K) > 4.5 appears to be 
better described by two power laws, rather than one, with 
the change of slope happening at 6 ~ 0.02 deg. The model 
also predicts such a change in slope, though this is milder 
than observed. The bottom panel in Fig. [TT] shows how this 
change in slope is due to the change from counting pairs of 
galaxies within the same halo, the one halo term, to counting 
those from different haloes, the two halo term, as was also 
seen for the spatial correlation function in U At scales larger 
than 9 > 0.2 deg we can see the change in the curvature of 
aj(6l) intrinsic to the two-halo term, that follows the shape 
of the dark matter clustering. 

Fig. I 111 shows that the model underpredicts the angular 
clustering with respect to observations. The difference be- 
comes larger at larger pair separations beyond 9 — 0.02 deg. 
However, as can be seen in Fig. 1111 the uncertainty due to 
sampling variance also increases with pair separation, where 
the clustering measured in the two fields shows the largest 
differences. 

As indicated before, the differences between model pre- 
dictions and observations are not due to the differences in 
the redshift distributions. Sampling variance could be be- 
hind the discrepancy between observations and model pre- 
dictions, at least at small scales, 9 < 0.02 deg. On larger 
scales the difference is clear even when allowing for the sam- 
pling variance errors in observations of different fields. On 
large scales, 9 > 0.02 deg, the angular clustering is under- 
predicted by a factor of 3 for EROs with K < 18.8. 

The match between predictions and observations is 
worse for EROs selected with (i — K) colour than it is for 
{R — K) colours. The photometric errors in the CTIO i-filter 
are larger than those for the r-filter, which implies that the 
selection made using {i — K) is slightly noisier than with 
(r — K), at least for the SA22 field. In fact, the prediction 
for EROs with (i — K) > 4.5 perfectly matches the angu- 
lar clustering estimated from observations of EROs selected 
with {i — K) > 4, potentially pointing towards a problem 
specific to the colour cut. The predicted angular clustering 
using the Subaru i-band is almost identical to that shown 
in Fig. 111! We have further explored the origin of this dis- 
crepancy by predicting the angular clustering of EROs with 
(i — K) > 5. We find a similar trend to that reported for the 
spatial clustering: the clustering of redder EROs is boosted, 
with the larger v ariation seen at small scal es. In line w i th this 
prediction, both lDaddi et all l|2000bl ) and lKim et ahl (|201ll ') 
observed that redder EROs have angular correlation func- 
tions with higher amplitudes. 

Optical colours in the model are affected by the treat- 
ment of the ram -pressure stripping in satellite galaxies 
jFont et al.ll2008l '). When selecting EROs by their {i ~ K) 
or {R — K) colours we could be enhancing their numbers by 
selecting satellite galaxies that are redder than they should 
be due to an oversimplification of their gas physics in the 



I Bower et al.l model. This would change the amplitude of the 
angular clustering, particularly at small scales. We have in- 
vestigated this p oint by calculatin g the angular clustering of 
EROs using the iFont et al. (12008 ) model. This model is an 
extension to that of iBower et al.l which introduces a more 
realistic treatment o f the stripping of hot gas in satellite 
galaxies. iFont et al.l also double the stellar heavy element 
yield, improving the match of the locus of the red sequence 
to observation al data at z ~ 0. 

Using the iFont et all (|2008f ) model we find that the 
number counts of EROs with (i — K) > 4.5 are overpre- 
dicted for A' < 18 by around a factor of two. The r edshift 
distribution of EROs predicted by the iFont et al.l model 
has a non negligible fraction of galaxies at z < 1, with a 
median redshift of 1.13, below that observed. At a given 
redshift, the spatial correlation function predicted by the 



Font et al.l m odel is very close to that predicted from the 
Bower et al.l model, for EROs with either {i — K) > 4.5 



or {R — K) > 5. Howev er, the angular clustering calcu- 
lated from the iFont et ahl model is predicted to be boosted 
by at most a factor of ~ 1.5 at lar ge scales, when com- 
pared with that from the lBower et al.l model. Therefore, the 
iFont et al.l model does not reconc ile the predi ctions with the 
observational estimates. Since the lFont et al.l model includes 
changes in both the stellar yield and the treatment of hot 
gas stripping th at have oppo site effects on galaxy colours, 
we have run the iBower et al.l model changing only the stel- 
lar yield to double the original value. T he results o f doing 
this exercise are similar to using the fuU iFont et al.l model. 
Both the redshift distribution and number counts change 
upon changing the yield. These changes are such that the 
angular clustering remains close to that obtained with the 
original model. Thus, we find that a change in stellar yield 
cannot account for the difference between the model predic- 
tions and the observations, and leads to some predictions 
actually agreeing worse with the observations. 



5.2.3 "Old" and "dusty" EROs with {i - K) > 4.5 

The number counts for EROs separated into "old" and 
"dusty" galaxies are shown in Fig. [TO] This separation has 
been done according to the loca tion of the EROs in the 
colour-colour diagram proposed bv lFang et al. I (120091 ). where 
the line ( J - K) = 0.20(i - A') -I- 1.08 is defined to set the 
boundary between the se typ e s. Th is selection is similar to 
that of .Pozz etti _& Ma nnuccil l|2000l ). but is tuned to the fil- 
ters u sed to select EROs with (i — K) colours. iFang et al.l 
([2009|) compared EROs selected with both this colour-colour 
selection criterion and on the basis of their spectra, finding 
that the two meth ods a gree reasonably well. Nevertheless, 
according to Fang et al.| . ~ 33% of the "old" EROs selected 
by the colour-colour method are actually young EROs when 
looking to their spectra. 

We find that about 90% of the EROs with A' < 19 are 
"old" galaxies according to the iFang et al.l colour criterion. 
We have tested that a small shift of about 6% in the bound- 
ary proposed by l^n^eFaLl, (J-K) = 0.2Q{i-K) + 1.02, 
produces a similar change in the percentage of "old" EROs 
of a 5%. Thus, the colour-colour separation is reasonably ro- 
bust for the predicted EROs, since they do not preferentially 
pop ulate the region c lose the colour-colour boundary. 

iKim et all (120111 ) observed that 60% of their EROs are 
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Figure 12. The predicted two-point angular correlation function 
of "old" (red solid line) and "dusty" (blue dashed line) EROs 
w ith (i — K) > 4.5 and K = 18.8. The observational estimates 
of lKim et all ll201ll) for "old" (circles) and "dusty" (stars) EROs 
are shown with Poisson error bars. See text for the details of the 
division between "old" and "dusty" EROs. 



"old" galaxies, a lower fraction than the mo del predi cts. 
This result is similar to the conclusion from IPaper j on 
the nature of EROs selected using {R — K) colours, al- 
though, in the case of the (i — K) selection the difference 
with observations is larger. One possible reason for this dis- 
crepancy is that although GALFORM broadly reproduces the 
observ ed colour bimodality observed for SDSS galaxies at 
z = l|Gonzalez et al.ll2009l ). the model predictions are dif- 
ferent in detail from the obser vations. However , we have 
fou nd that when using either the lFont" et all (j2008l ) model or 
the iBower et al.l model with twice the original yield value, 
about 60% or 67%, respectively, of the predicted EROs with 
K ~ 18 are "old", in good agreement with observations. 
Therefore, the difference in the split between "dusty" and 
"old" ERO s appears to be connected with the value of the 
yield in the iBower et al.l model. 

Fig. [To] shows that the percentage of "old" EROs starts 
to drop for fainte r magnitude bins . The same tendenc y has 
been observed bv lKim et aD (|201l[ ). iKong et all (j2009h also 
found the same trend for a sample of EROs selected with 
(7 - ifs) > 4 and Ks < 19.2. 

The predicted median redshifts f or "old" and "dii sty" 
EROs are 1.44 and 1.25, respectively. iKim et al.l (j201ll ) es- 
timated the median redshifts for the observed "old" and 
"dusty" EROs to be z = 1.36 (10% and 90% percentiles: 
z = 1.1 and z = 1.7, respectively) and z = 1.25 (10% and 
90% percentiles: 2 = 1 and z = 1.8, respectively), respec- 
tively, in reasonable agre ement with the m odel predictions. 
Contr ary to these results. fPang et al.l (|2009l ) and Kong et"al] 
(|2009l ) found more "dusty" EROs at higher redshifts. How- 
ever, we have found that this statement is strongly depen- 
dent on the particular {i — K) cut made, since when using 
(i — K) > 5 we find that "dusty" EROs are predicted to 
appear at higher redshifts than "old" ones. 

Fig. [12] shows the predicted two-point angular corre- 
lation function of EROs with {i — K) > 4.5 separated into 
"old" and "dusty" galaxies. Fig.[T2]shows that the predicted 
clustering for both types of EROs differ only on small scales. 



Thus, model predictions and observations agree in that both 
"dusty" and "old" EROs are clustered very similarly. 

In the previous section we discussed the possible ori- 
gin of the difference between the predicted and the observed 
angular clustering for EROs with (i — K) > 4.5. The bal- 
ance between "dusty" and "old" EROs when using {i — K) 
colours could alter the predicted angular clustering. How- 
ever, Fig. [12] shows that the clustering of "dusty" EROs is 
both observed and predicted to be similar, though slightly 
less strong than that for the "old" EROs. Therefore, a larger 
fraction of "dusty" EROs within the model would actu- 
ally slightly increase the difi'erence between observations and 
model predictions, by further reducing the predicted angular 
clustering at large scales. 



6 CONCLUSIONS 

In this paper we have extended the tests of the 
GALFORM galaxy formation model by continuing the 
study of Extremely Red Objects (EROs) started with 
iGonzalez-Perez et all (|2009D (Paper I). EROs are massive, 
red galaxies at 0.7 ^ z < 3 and their numbers and proper- 
ties have posed a challenge to hierarchical galaxy formation 
models. In this paper we have analysed the halo occupation 
distribution (HOD) and clustering of EROs predicte d based 
on the published model of galaxy formation of Bow er et al.l 
(|2006l ). The parameters in this model were set to repro- 
duce observations of the local galaxy population. No pa- 
rameter has been changed for the work presented here. The 
iBower et al.l model gives an impressively close match to the 
number counts of EROs (Paper I). This model also matches 
the evolution of the K-band luminosity function and the in- 
ferred evolution of the stellar ma ss function. 

The HOD predicted by the H ower et al] model has a 
strikingly different form from the canonical assumption of a 
step function that reaches unity for the central galaxies plus 
a power law for satellites. The central galaxy HOD has a 
somewhat more rounded form than a step function and does 
not reach unity for K-selected galaxies or EROs. This is due 
to the non-monotonic relation between host h alo mass and 
central galaxies luminosity. According to the iBower et al.l 
model, at 2 = 2.1, galaxies with K < 20 are very bright, 
L > L*. For such central galaxies, the AGN feedback in this 
model truncates the initially monotonic mass-luminosity re- 
lation and increases the scatter, modifying the shape of their 
HOD. Our predictions suggest that a revision is needed to 
the canonical form assumed for the HOD to model the clus- 
tering of bright galaxies. 

The HOD of EROs suggests that the {R - K) or {i - 
K) colour cuts do not select all galaxies above a certain 
mass threshold since they leave out some of the younger 
star forming galaxies. However, EROs with {R — K) > 5 
or (i — K) > 4.5 come very close to representing the whole 
population of galaxies with A" < 20 at 2 > 1.5. On average, 
halos more massive than 10^^h~^ Mq host at least one ERO 
with K <2Q and (R - K) > 5. 

At a given redshift, the contribution of satellites dom- 
inates the shape of the HOD of the redder EROs. The 
satellite contribution also becomes incr easingly domi nant at 
lower redshifts. We predict that, in the lBower et al.l model. 
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brighter and redder EROs should be found in more massive 
haloes. 

We have found that the minimum halo mass needed to 
host two EROs can be more than an order of magnitude 
larger than that needed to host just one ERO. A similar 
factor is found for K-selected galaxies at 1 < 2 < 2.5, so 
this result is not peculiar to EROs but rather is intrinsic to 
bright galaxies. 

We predict from the lBower et al.l model that the spatial 
two-point correlation function of EROs is roughly consistent 
with a power law for scales 0.03 '^^t <iQ ^~^Mpc, though in 
real-space there is a ~ 10% change of slope at r ~ 1 /i~^Mpc, 
which reflects the shift in dominance from the one-halo to 
the two-halo term. EROs with (7? - isT) > 5 and K < 20.2 
at z = 2.1 are predicted to have a bias of ~ 3. 

The predicted clustering does not depend strongly on 
apparent magnitude over the range 18.4 < K < 20.2. We 
predict the clustering of redder EROs to be boosted, partic- 
ularly on small scales, due to a higher proportion of satel- 
lite galaxies. We have found that those EROs predicted 
to be experiencing a star formation burst are more clus- 
tered than those which are passively evolving. However, 
when separating EROs by their specific star formation rate, 
SFR/M = lO^^^yr^^, active EROs are predicted to have a 
correlation length of ro = 5.7/i.^^Mpc, which is lower than 
that predicted for quiescent ones, ro = 7.7/!.~^Mpc. The lat- 
ter value is comparable but lower than the observational 
estimations. 

The predicted spatial correlation length of EROs is 
smaller than most observational estimates. However, since 
the match between the predicted angular clustering and the 
observed one for EROs with (7? — K) > 5 is rather good, it 
is quite likely that the discrepancy in the correlation length 
is due to the assumptions required to derive this quantity 
from the observations. Most of the observational studies of 
EROs rely on photometric redshifts, which introduce a large 
uncertainty into the redshift distribution, which is needed 
to obtain the spatial correlation function from the angular 
one. The surveyed area usually is quite small and the inte- 
gral constraint has to be included to account for the bias in 
the inferred mean galaxy density and the impact this has on 
the estimated clustering. Another simplification that is usu- 
ally made is to assume that the spatial correlation function 
is a pure power law. Our model predicts departures from a 
power law, in agreement with observations of galaxies at low 
redshifts. Often it is also assumed that the clustering varies 
monotonically with redshift, and this evolution is usually 
parametrised as a global, scale independent factor. However, 
our results suggest that the evolution of ^ with redshift is 
not monotonic and furthermore that it varies with pair sep- 
aration. 

The predicted angular correlation function for EROs 
with {R—K) > 5 matches the observational estimates within 
the range of error allowed by sampling variance. On a scale 
of ~ 0.02 deg, we find a small change of slope in the 
angular clustering due to the transition from the one-halo 
to the two-halo terms. 

We have also explored the predicted angular clustering 
of EROs selected by their (i — K) colour. This colour se- 
lection is more efficient at leaving out galaxies with z < 1. 
We find that, unlike the case of EROs with (R — K) > 5, 
the angular clustering of EROs with (i — K) > 4.5 is under- 



predicted, particularly at large separations. We have shown 
that this difference is due partly to sampling variance. How- 
ever, at large scales sampling variance alone cannot explain 
the factor of 3 difference between the model and the obser- 
vations. We have explored the possible origin of this discrep- 
ancy. We have shown that the impact of using the redshift 
distribution estimated directly from observations, using dif- 
fere nt i-band filt ers or doubling the default stellar yield in 
the I Bower et al.l model has a minimal impact on the pre- 
dic ted angular correlation at large scales. We have also run 
the lFont et al.l ("2008) model, finding that, on large scales, it 
predicts an angular correlation f unction only slightly higher 
in amplitude than that from the lBower et al.l model. Unlike 
in the observations, EROs with {i — K) > 4.5 are predicted 
to be dominated by "old" galaxies, classified as s uch by their 
location in the colour-colour space proposed bv iFang et al.l 
1 20091 ). However both the predicted and the observed clus- 
tering for "old" and "dusty" EROs are very similar, and 
thus a change in the split between these two types cannot 
account for the discrepancy in the angular correlation func- 
tion on large scales between model and observations. Nev- 
ertheless, the problem appears to be related to the i-band 
filter, since our prediction for EROs with (z — K) > 4.5 per- 
fectly matches the observed angular clustering of EROs with 
{i-K)> 4. 

Overall, the predictions based on the lBower et al.l model 
match observations reasonably well, once sample variance in 
the current data is taken into account. 

This is the second paper in a series which examines the 
properties and nature of red galaxies in hierarchical models. 
In the third paper we will examine different colour cuts used 
to select red galaxies at 2 > 1 and compare the properties 
of their present-day descendants. 
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